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Based on bond arguments, a Hamiltonian is introduced to describe the fundamental 
physics of the collective rotations of oxygen atoms in oxides. Values for the relevant 
material parameters are estimated for silica and Cu oxides. Zero- temperature 
phase diagrams are presented; the different phases and associated transitions are 
characterized. Phases with non-phononic excitations are found. 



1 Introduction 

The tendency of oxygen atoms to form just two bonds with a bond angle 
smaher than 180° allows the appearance of low energy atomic motions related 
to internal rotational degrees of freedom which, in a first approximation, in- 
volve neither bending nor stretching of bonds ("rigid unit modes" in silicates, 
"floppy modes" in glasses) □ A huge variety of compounds exhibit the oxygen- 
bridge X-O-X basic unit, X being a linking atom. Among the many silicates, 
/3-cristobalite (a silica polymorph) shows these internal rotations most clearlyo 
The perovskites are also interesting candidate materials, especially the Cu02 
planes of the high Tc superconducting materials. 

2 Model and Hamiltonian 

A lattice of fixed X atoms is considered with O atoms bridging between them. 
The degrees of freedom are the angles 4>i associated with the oxygens in their 
rotations around the axes defined by the neighboring nodes (given a fixed X-0 
distance those rotations are the only possible motions). The Hamiltonian is 

n = -B^d^/dcfyl + 2t ^ cos{(j), - - 2C^ cos(n0.O> (1) 

including the kinetic energy, a nearest neighbor interaction which tends to 
align the O bridges {t > 0), and a hindering potential. The latter accounts for 
effects in real solids which impede the rotation of the oxygen atoms. It is a 
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Figure 1: (a) Dependence of t and B on the puckering angle a, (b) The energy barrier Ei, 
and lowest vibrational energy vs. a for the radial motion of oxygen. Si02 case. 

fixed potential, whereas in the real solids the impeding objects (other atoms) 
are mobile, n gives the number of potential minima around the circle. 

The rotor constant is B — h^/2I, where I — moR^, and i? = 6sinQ;, b being 
the bond length (1.6 A for Si-0, 1.9 for Cu-0), and the puckering angle a — 
(7r-2 9x-0-x)/'^- Therefore, B^B^/ sin^ a, with Bo « 0.06 meV for Si and 0.04 
for Cu. The interaction parameter t accounts for the energy needed to rotate 
an O atom while keeping its neighbors fixed. This energy is essentially used to 
bend the 0-X-O angle away from its natural value. Then, t = to{\ — cos 2a). 
Using the 0-X-O bond bending force constants, to « 170 meV for Si, and 
230 for Cu. B decreases with a whereas t increases, defining an a for which 
B =^ t (7° for Si, 5° for Cu). C is much more dependent on particular crystal 
structures and more difficult to estimate. It increases with a,|-|and can be as 
small as 0.1 meV for a — 20°, as found for interstitial O in GeB 

The radial vibration of the O atoms involves the stretching of X-0 bonds 
and can be neglected for well-defined rotor situations. However, for small 
enough a, the radial motion becomes important since AR ^ (R). An estimate 
of the range of validity of H can be drawn from the X-0 bond stretching force 
constants, by comparing the ground-state energy for the radial vibrations 
with Eb, the potential energy barrier at i? = 0. Ef, = K{1 — cosa)^, with 
« 80 eV for Si and 30 for Cu. E^ = K' sin a, with K' « 70 meV for Si, 35 
for Cu. Eb > El for a > 9° both for Si and Cu. Only the (fa will be coi^idered 
here, the effects of the radial motion being the subject of future worko 

3 Phase diagram 

For C = we have a one parameter {t = t/B) problem. The zero-temperature d- 
dimensional quantum Hamiltonian lies in the universality class of the classical 
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Figure 2: Dispersion relation for the elementary excitations for the ID Hamiltonian, (a) for 
t<^B, and (b) for t^B. (c) shows the dependence of the gap A on t. W, a, and Vs stand 
for band width, lattice parameter, and sound velocity, respectively. 




{d+ l)-dimensional XY modela This model exhibits a phase transition at some 
critical tc of the order of 1. It is an order-disorder transition for d> 1, implying 
a symmetry breaking in </) for the ordered phase {t > tc), with its associated 
long-range order and gapless Goldstone excitations. For d=l the transition 
is of the Kosterlitz-Thouless kind. The i > ic "quasi-ordered" phase displays 
algebraic order ((e*^"^'"'^^)) ^ — with ij < 1/4). The disordered phase 
shows an exponential decay of (e'('^'~"^j)), with gapped excitations. For d=l, 
near the transition, the gap closes like A '--^ 1/^ ~ ^-K/y/ia-i^ ^ being the 
correlation length [Fig. 2(c)]. 

The limits help to understand the nature of the phases. For t<^l, the basis 
of independent-rotor eigenstates provides a natural language. The interaction 
term acts like a raising-loweriiig operator on the rotor quantum numbers m^. 
The original Ti can be mappecfl into a Hamiltonian of bosonic particles and an- 
tiparticles (m positive or negative) with a chemical potential fi = B, a, Hubbard- 
like repulsion U = B, and a hopping energy t. Using a Bogoliubov transforma- 
tion, the t limit can be solved yielding an independent-rotor-like ground 
state with quantum fluctuations of particle-antiparticle pairs, and rotational 
excitations (rotatons, dressed rotor excitations), with = Byl — Aicoska 
[Fig. 2(a)]. For t3> 1, a harmonic approximation of the interaction term yields 
gapless spin- wave-like excitations [Fig. 2(b)]. 

Information for the full phase diagramJ^C 7^ 0) can again be obtained 
from the corresponding classical d+1 modelQ For d=l the phase diagram is 
displayed in Fig. 3 (a)-(c). In addition to the phases discussed above, a third 
phase appears, an ordered phase corresponding to the locking of around one 
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Figure 3: Phase diagrams in the C—t plane for different values of n, the number of minima of 
the hindering potential, (a), (b), and (c) are for one dimension, (d) is obtained from a mean- 
field approximation, dependent on the coordination 2; the phase boundaries for n=2,4,6 are 
indicated in the graph. In every diagram the phase at the left is a disordered phase; the 
ones at the right, up and down, are ordered phases; stars denote critical points. 



of the n equivalent minima. For n > 4 the critical region separating the ordered 
and disordered regions has a finite width; it shrinks to a line for n = 4. For 
n = 2 the model belongs to the Ising universality class (C 7^ 0). Indeed, the 
Hamiltonian maps into a Ising model in a transverse magnetic field for C^B. 

A mean-field calculation has been performed to obtain quantitative esti- 
mates. We approximate the interaction term by 4zt(cogLiA)^^ cos z being 
the coordination number and (cos (j>) the order parameter^ The phase diagram 
[Fig. 3(d)] is expected to be more accurate in higher dimensions. The critical 
regions disappear, with classical order-disorder phase transitions remaining. 
Higher n increases the C range of the disordered phase. 
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